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Abstract. We study the problem of recovery the source a{t,x)F{x) in the wave equation in anisotropic 
medium with a known so that a(0, x) ^ with a single measurement. We use Carleman estimates combined 
with geometric arguments and give sharp conditions for uniqueness. We also study the non-linear problem 
of recovery the sound speed in the equation utt — c^{x)Au = with one measurement. We give sharp 
conditions for stability, as well. An application to thermoacoustic tomography is also presented. 



1. Introduction 

The purpose of this paper is to give sharp conditions for a recovery of a source term in the wave equation 
in anisotropic media modeled by a Riemannian metric by a single boundary measurement. In the process we 
give a more geometric treatment of the problem. This linear problem appears as a linearization and actually, 
as the full non-linear version, of the problem of recovery a sound speed, given the source. It has applications 
to thermoacoustic tomography. We are also inspired by the related works flj [5l |6j , [H Theorem 8.2.2]. The 
method in these papers uses Carleman estimates for hyperbolic inverse problems that originates in the work 
of [5] who considered the case of the wave equation with a potential with non-zero initial data. 

The main problem that we have in mind is the following. Let c{x) > and let u solve 

( {d^-c^A)u = in(0,r)xR", 

(1) < ^^|t=o = /, 

[ dtu\t^o = 0, 

where c = c{x) > and T > is fixed. Let c = 1 outside some domain fl with a smooth strictly convex 
boundary. Given /, and u restricted to [0,T] x dH., is it possible to reconstruct the speed c? Ideally, we 
want to do this with data on a part of dfl, as well. Next, assuming that we can, how stable is this? Clearly, 
some conditions on / are needed since when / = 0, for example, we get no information about c. This inverse 
problem is clearly non-linear. 

If we have two speeds c and c, then w ^ u — u solves 

( {df-c^A)w = a{t,x)F{x) in(0,T)xR", 

(2) I w\t^o = 0, 

[ dtw\t=o = 0, 

with 

(3) F:^d'^~c^, a^Au. 

We consider the more general linear problem of recovery of a function F, given a and w restricted to 
[0, T] X 9f2, or on a part of it. Again, some condition on a is needed since when a = for example, F cannot 
be recovered. We actually replace c^A in ^ by the Laplace-Beltrami operator Ag related to some metric, 
plus lower order terms. 
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Similar problems but for the recovery of a potential or the term p{x) in V • pV are studied in [H [51 [B]. 
A general abstract theorem of this type is presented in jB] Theorem 8.2.2]. In [B] and [5j, the principal 
part of the wave equation has variable coefficients, thus the geometry is non-Euclidean. This requires some 
assumptions on the speed or the metric. The method of the proof is to use Carleman estimates, and the 
assumptions are needed to satisfy the pseudo-convexity condition. Those assumptions are not sharp however. 
In fact, the proofs are essentially "Euclidean" , and roughly speaking, the conditions on the speed or on the 
metric require that the proof still works under an Euclidean treatment. Also, one global pseudo-convex 
function is used. One of the goals of this work is to sharpen those conditions thus extending the results to 
a larger class of speeds/metrics, formulate them in a geometric way and also prove local results. 

There are many works on related problems, including boundary control, for example, [U [TOl [22l |2T] . The 
conditions on the metric there are more geometric, requiring existence of a global convex function, or a 
somewhat general condition of existence of a global vector filed with certain properties. The proofs are 
still based on Carleman estimates but the goal is to recover non-trivial initial conditions, assuming, say, 
a Neumann boundary condition, and measuring Dirichlet data on a part of dfl. The conditions on T are 
formulated in terms of lower bounds of the speed and are not sharp. The geometry of the rays in those 
problems however is different than the application that we have in mind — there are reflections at the 
boundary. On the other hand, the methods there could probably be adapted to the problems studied in the 
works that we cited above. 

The pseudo-convexity condition needed for the Carleman estimates that we use is satisfied by assuming 
that the region where we prove unique continuation is foliated by a continuous family Eg of strictly convex 
surfaces. In case of data on a part F of dft, we require those surfaces to intersect dft in F. In contrast to the 
other works in this direction, we are not trying to construct one strongly pseudo-convex function. Instead, we 
prove unique continuation by incremental steps, each time using a different strongly pseudo-convex function. 

We describe the results in the paper now. We start in section [2] with the uniqueness Theorem 2.1 that 
is a version of 13 Theorem 8.2.2]. The time interval is (— T, T), there is no initial condition for wt a.t t = 0, 
and we study the problem of unique recovery of in ([2| given Cauchy data on a part of {—T, T) x dVl. We 
view this theorem more a as a tool than a goal, and the requirement that the time interval is symmetric 
about t = will be satisfied later by studying problems with solutions that have even extensions in t, like 
Q. In the rest of that section we show two ways to satisfy the convexity requirement. First, assuming that 
dVl is strictly convex, we show in Theorem 2.2 that F = in some collar neighborhood of dO, of the kind 
dist(a;, 9r2) < T ^ 1, as long as the surfaces dist(x, 957) — s, s € [0,T] are smooth and still strictly convex. 
The second one is to show that F = in a subset of f2 that can be foliated by strictly convex surfaces starting 
from ones outside see Theorem |2.3[ This only requires Cauchy data on a part of 9f2, where those surfaces 
intersect 9r2, and proves F = in a subdomain. The condition on T is sharp. We give a few examples. 

In sectionjsjwe study the non-linear problem of recovery of the speed c in ([T]) and the linear one of recovery 
of the source F in ([2]) described above. The time interval is (0, T) now, but the initial condition = for 
f = in ([I]), and the requirement that a in ([2]) has a sufficiently regular even extension in the t variable 
allow us to use the results in the previous section. In contrast to the boundary control problems, in the 
thermoacoustic problem we are given the Dirichlet data on (0,r) x dU, or on a part of it but no Neumann 
data. On the other hand, we know that the solution extends for a; ^ f2 as a solution again, and the initial 
data at t = is zero there. This allows us in Lemma [3.11 to recover the Neumann data from the Dirichlet 
one in case of data on the whole dO.. Then we extend the solution in an even way for t < Q and apply the 
results in section [2] The main requirement is 57 to have a foliation of strictly convex surfaces, and the time 
interval (0, T) is sharp. 

The partial data case, with observations on (0, T) x F, where F C 957 in the thermoacoustic setting is 



studied in Theorem 3.2 The main difficulty is the need to recover the Neumann data there as well. Then one 
applies directly the results in section[2j We show that one can recover F, respectively c in some neighborhood 
of F that might be smaller compared with the case of having Cauchy data on the whole (0, T) x F. There is 
a new "cone" condition that might shrink the domain where we prove F = 0, or respectively c = c. 

At the end of section |3j we study the stability of the linear and non-linear problems for ([T]) and Q , 
respectively. As a general principle, for stability, we need to be able to detect all singularities, see ((5l|). For 
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the Hnear problem at least, this is a necessary and sufficient condition for stability in any Sobolev spaces, 
see |13j . The corresponding stability estimates are formulated in Theorem 3.4 and Theorem 3.5 
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2. A UNIQUENESS RESULT FOR RECOVERY OF A SOURCE WITH ONE MEASUREMENT 

Let f2 be a bounded domain in R" with a smooth boundary, and let 

(4) P = 52-A, + ^6,9,.+c 

be a differential operator in Q :— (~T,T) x C R", where g is a smooth Riemannian metric on Q, and Uj, 
b, c are smooth functions on Q. 

The level surface E = {V' = 0} of some smooth function ip is called strongly pseudoconvex w.r.t. the 
hyperbolic operator P, if 

(i) E is non-characteristic, i.e., {tptl \d.x''P\ when "0 = 0, and 

(ii) H^ip > on r*f^ \ whenever ip = p = Hpip = 0, 

where Hp is the Hamiltonian vector field of the principal symbol p = — + of P, and A is the dual 
variable to t, see, e.g., [T5]. Here and below, | • | is the norm in the metric g of a covector or a vector. The 
second condition says that E is strictly convex w.r.t. to the null bicharacteristics of p, when viewed from 
ij} > 0- In other words, the tangent null bicharacteristics to E are curved towards ip > 0. 

The function (j) with a non-degenerate differential and non-characteristic level set = is pseudoconvex, 
if a condition stronger than (ii) is satisfied. We are not going to formulate that condition; it would be enough 
for our purposes to use the well known fact that if the E = {t/; = 0} as above is pseudoconvex, then for 
fi ^ 1, (f) = exp{fiip) — 1 is a pseudoconvex function, non-degenerate on E, and E = {0 = 0}; moreover 
{(t> > 0} = {ip > 0}. For details we refer to [T^ . 

Let (j> be strongly pseudoconvex in Q w.r.t. P. Then it is well known that for all u e Cq{Q), 

(5) T J e^^''' {u^t +\'^u\'^ + T^\u\^) dtdx <C J e^^'^\Pu\^dtdx, r > 1, 
see [HIIHIIZ]. 

To reformulate condition (ii) in the tangent bundle, recall that the metric g provides a natural isomorphism 
between covectors and vectors by the formula <& : {x,£_) >->■ {x,v), where = gij{x)v^, in local coordinates, 
where w is a vector at x. For any function ijj on T*n, one gets a function on Tfi. Let q = |'CP/2 be 
the "x part" of p, rescaled for convenience. It is known that Hq = <1>*G<I>*, where G is the generator of the 
geodesic flow. Also, the energy level q — 1/2 is pushed forward to the unit sphere bundle Sil. 

We have Hp/2 = —Xdt + Hq, therefore, 

Identify covectors with vectors by the map $, to get that condition (ii) is equivalent to 

(6) 1^ = 0, V't - GV^ = =^ {dt ~ Gfi, > for (t, x) G Q, |^| = 1, 

and we used the fact that p — implies A^ — as well as the homogeneity properties of G w.r.t. ^. 
Let us look for -0 of the type 

^ = r'^{x)- 5t^ - s, 0<(5<1, 
with s a parameter, Then to satisfy (ii), it is enough to have 

(7) G\r^/2) > S\e- 
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Since we want eventually to take the limit (5 — > 1 to get sharp results, we arrive at the condition 

(8) G^(rV2) > 

It is enough to have this inequality in the a;-projection of S in f2, as the first condition in ([g]) indicates. Note 
that this guarantees (ii) only, we still have to choose r so that (i) holds. The latter is equivalent to 

(9) \A{r^/2)\^5\t\ on Eng. 

Example 1. Let r{x) — \x — xq\ in R" with some fixed xq] then ijj = \x — — 5t^ — s. Then 

G\r^l2)^{^-V,f\x~xo\^l2 = \^\\ 

and Q holds. Condition ^ is satisfied for s > because \d{r'^/2)\^ ^ \x - xq\^ = 5t^ + s > . Then S 
is a hyperboloid of one sheet. 

More generally, let r{x) = p{x,Xo), where p is the distance in the metric g. This is the function that 
has been used in the Riemannian case. It satisfies Q for r ^ 1 only, in general. For metrics of negative 
curvature, there is no restriction, see [T?l . 

In this paper, d/dv denotes exterior normal derivative to dU,. 

Theorem 2.1. Let Q = {-T,T) x fi, and 

(10) dia e C{Q), j<2, Fe L^{n), and dld"u e L^{Q), j < 3, \a\ < 1. 
Let u be a (non-unique) solution of 



(11) 



Pu = a{t,x)F{x) in{~T,T)xn, 
u\t=Q = in n 



Let (j) he a strongly pseudoconvex function in Q. Let 

(12) g C {-T,T) X dn, (l)<0 on {{-T,T) X dn)\g, and(j){t,-)<(j3{0,-) for\t\<T. 
Let supp F C K , where I'C d Q is compact, and 

(13) a(0,-)^0 onK. 
If 

(14) u = du/dv = ong, 
then 

F = in {x; (j>{0,x) > 0}. 

Proof We follow the proof of [H Theorem 8.2.2]. Set Q^ = Qr\{(f>> e}. Fix e > 0, and let x e C°° be such 
that X = 1 in Q^, and suppx C Qq. We will apply the Carleman estimate (jsj to dfxu, j = and j — 2, by 
shrinking Qo to on the left. We are using the fact here that u has zero Cauchy data on g. Then 9jX" 
can be approximated by C^{Q) functions in the norms, see also the remark following the proof. We 
have 

Pdlxu = dl {xPu + [P, x]u) = dl ixaF + [P, x]u) , 
where the commutator [P, x\ is a differential operator of order 1. Since x 1 on Qc, we get 



(15) T 



/ e'^^{rM' + r'\uu\' + \u,unda<c( [ e'^^\F\' da + [ e^^^ ^ |a,^9,>p da 



< C I e2"*|P|2da + Ce2"% 
Q 
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where da = dtdVol(a;) is the natural measure on Q. We will estimate the first term in the r.h.s. above. 
From equation (11) and is initial condition, utt(0, •) — a(0, •)/. By the assumptions on a, < C\utt{0^ ■)\. 
By ([12]), 

(16) / e^^^lFpda < 2r / e'"'^(°'- Vtt(0, •)!' d Vol(x). 

Jq Jn 

This integral admits the estimate 

s2-'^(°')|u«(0,.)pdVol(x) = - r / I- 

Jo Jn OS 



(e2-^(^--)|w„(s,.)P)dVol(x)ds+ / e2-^(^'-)|u«(T,.)pdVol(a;) 



<C [ e^-^ (r|Mtt|2+T-V„t|2)da + C, 
Jo 



because (/'(T, ■) < 0, and by the Cauchy inequality. This inequality, together with (16), estimate the integral 
of the first term in the r.h.s. of (15). Therefore, 



• / ^r^{Au? + Auu\'' + \uuA'') da<cfr-i / {Autt\^ + \uut?) da + , 



,2te 



Split the integration on the right into and Q\Qe- For r ^ 1, the integral over Qe will be absorbed by 
the l.h.s. On Q \ Qe, we have e^^"^ < e^^^ Therefore, 

r / e^^-^lwl^ da < Ce^^", for r > 1. 

Jq. 

Thus, 

/ e^^^-^-") |wp da < C/t, for r > 1. 

Jq. 

Since cj) — e > in Q^, letting r 
proves the theorem. 



oo, we get M = in Qe. Since e > is arbitrary, we get m = in Q. This 

□ 



The following lemma will allow us below to apply the proof of the theorem to a larger class of non-smooth 
boundaries. 

Lemma 2.1. Let D C R" be open, and assume that near each point xq G dfl, D is represented by y > 
and z > 0, where y, and z are Junctions with non-zero differentials, satisfying the following condition 

if {y = 0} and {z — 0} intersect, then {y = 0} \ {z = 0} is dense in {y — 0}. 

Let u G C^{D) has Cauchy data on dD in the sense that extended as outside D, it still belongs to . 
Then u can be approximated by Cq{D) functions in i7^(R"). 

Proof. Let x £ ^'^"(R) be such that x = 1 near 0. Set locally, near xq £ dD, 

= (1 - X(2//£))(1 - Xiz/£))u. 

Using a partition of unity, we define such G C^(D). We claim that Me u in iJ^(R"), as e 0+. 
Take the second derivatives of u — We to see that we need to show that the following terms converge to in 
L2(R"): 

e~\'iy/e)u, e-\"(y/£)u, s-^x'iy/e)x'iz/e)u, 
where we used the fact that the derivatives of y and z are bounded. We list terms involving lower powers 
of and derivatives of u, which analysis is similar. Similarly, we will not analyze the zero and the first 
order derivatives of u^. Since y/e and z/e are bounded on the support of x'(y/e) and x'iz/e), respectively, 
we may replace the leading coefficient e^^ in the first term by y^^, etc. Since m = for z = 0, and dy ^ 0, 
we have u — zu\, where U\ is a smooth function. Next, ui = for y = at least when z ^ 0. That set of 
y's however is dense in \y = 0} by assumption. By continuity, wi = for ?/ = 0. Thus u\ — yu2 with U2 
smooth, therefore, u — zyu2- Now, the proof of the claim follows from the fact that x"{yl^) tends to in 
L^(R"); and this is also true if we replace x" by x' or x- D 
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We recall next a unique continuation result due to Tataru [TH], see also [T^ Theorem 4]. Assume that a 
locally function u solves the homogeneous wave equation Pu = (near the set indicated in (17) below) 
and vanishes in a neighborhood of [—T,T] x {xq} for some xq and T > 0. Then 

(17) u{t,x)=0 for |t| +dist(a;o,a;) < T. 

One can formulate a uniqueness continuation statement along a curve in the x space. 

Proposition 2.1. Let [0,T] 9 cr i— > c(s) be a smooth curve in R" so that c(0) = xq for some xq, and a is a 
unit speed parameter in the metric g. Let u — near [— T. T] x {xo}. Then 

(18) u{t,x) = near {{t,x); x — c{a),\t\ < T — a}, 
provided that u solves the homogeneous wave equation Pu — near that set. 

The proof of the proposition is contained in the proof or il6, Theorem 6.1]. The idea is to cover the curve 
c with small balls and to prove unique continuation step by step. 

Next theorem in fact follows from Theorem 
the proof of Theorem |2.3| We refer to Figure [T 



2.3| below but its proof is simpler, and serves as a basis for 
for an illustration. 



Theorem 2.2. Assume that dfl is strictly convex. Let T > be such that x" :— dist(x, dQ) is a smooth 
function in f2 with non-zero differential for < x"^ < T; and {x" — s}, < s < T are strictly convex 
surfaces. Let u solve (11) and the function a satisfies (13) for < a;" < T. Assume also the regularity 
conditions (10). Then if u has zero Cauchy data on {—T,T) x dQ, we also have 



(19) 



F{x) =0 for X en, dist(2:, dn) < T. 



Proof. Let {x',x"') be normal boundary coordinates near dfl with sq fixed. Here x" is the signed distance 
to dil so that x" > in il. The function x" is defined in a small neighborhood of dil while x' are local 
coordinates near some boundary point. The metric g then takes the form 

5„^(a;',a;")dx"da;'^ + (da;")2, 

a, /3 < n — 1. Given f{x), independent of the dual variable we have 



(20) 



G'f=U 



dx^ 



. d 



e—f 



d'f 

dx^dx^ 



a/ 

dx'- 



Here, r"^, — —^dgjk/dx"- is the second fundamental form H > of the level sets of x", w.r.t. the chosen 
orientation, and it is zero when either i = n oi j = n. Assume now that / is a function of x" only. Restrict 
^ todn = {x" = 0} to get 



d'f 

a(x"): 



(r)'-n(x')(r,r) 



9/ 

9x" 
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To satisfy ^ for / 
(21) 



we need 



> R for x''' 



0, 



d^f df 
9(x")2 - 'dx^ 

where 1 /R is the minimum over dfl of the smahest eigenvalue (principal curvature) of the second fundamental 
form II. We can think of R as the largest curvature radius of dfl; and by assumption, < i? < oo. Then 
the following function satisfies ([s]): 

r{x) = 

because then f :— = {R~ 



) /2 clearly satisfies (21). Therefore, the function 



(22) 



,ri\2 



Mt,x) := (i?-x") 



St' 



is strongly pseudoconvex, assuming also that = is non-characteristic. Also, the last inequality in ( 12 ) 
holds. Note that in Example [T] if dfl — {x; \x\ = Rq}, then a;" — Rq — \x\, and one can choose R — Rq- 
Then ijj = — 5t^ — s, which is the phase function in the example. As in the example, we can show that 
{ips = 0} is non-characteristic for s > 0. In fact, we will work locally near a:" = t = 0, and s = R, and 
clearly, ips is non-characteristic there. 

Fix < e <C 1, and let < < e. We restrict s to the interval (i? — e)^ < s < i?^. This choice of the 
phase function corresponds to pseudoconvex surfaces given by 

(23) {R-x^'f -St^ ^s, [R-ef <s<R^, < x"" < e < R. 

In Rt X Clx, this restricts t to \t\ = 0{y/e), see also (24) below. We will apply Theorem 2.1 with a the phase 
:= exp(^?As)-l, ^ > 1, with s as in (23). On g := {-T,T) x dfl, we have V'sU"=o = R^-St'-s. 
< outside Q, we need -0^ < there, and therefore R^ — ST' < s for all s as in (23). Therefore, 



T> ^R^ -{R- £)2 = 



function ( 

To have ( 
if 

(24) 

always true when e ^ 1, we can apply the theorem for 6 < 1 close enough to 1. Therefore, we get that if u 
has zero Cauchy data on (— T, T) x dV. with T as in (24), then F{x) — for "08(0, x) > for any s as before, 
i.e., for < x" < £. Note that this does not prove the theorem yet, even when T is small enough so that we 
can have equality in (24) with some e satisfying the smallness requirements. The reason is that we get F = 
in a much smaller region: < x" < e instead of < x" < T because for small e, we have T ^ %/e ^ e. 
Also, r ^ 1 when R is large, i.e., when 917 is close to a flat surface at some point and direction. In other 
words, the price that we pay with T to push suppf by e depends on the (largest) radius of the curvature 
of 9ri, and this is not what we are trying to prove. We will use this argument as an incremental step only, 
and will prove the theorem by applying a finite number of such steps. 

To get the sharp time T, not necessarily small, we notice that we just proved that if u has zero Cauchy 
data on (-T, T) x dVL with T as in ([24| and e < 1, then F{x) = for a; € f), dist(a;, 917) < e. Then Pu = 
in the same domain and \t\ < T, by (11). By unique continuation, see Proposition 2.1 

(25) u{t, x)^0 for a; e n, dist(a;, dn) + \t\ < T, dist(x, dn) < e. 
In particular, 

(26) d^uit,x) = for X en, dist{x,dn) ^ s, \t\ <T - e, \a\<l, 
provided that T > e. 

Let £ be the supremum of all e for each the following statement holds: if u has zero Cauchy data on 
{—T,T) X dfl, then F{x) — when dist(a;, 912) < e. Then e has that property as well. If £ < T, then by 
the argument above, see ([26]), u has zero Cauchy data on (— T + £, T — e) x 912. Then we can repeat the 



argument leading to (25) to reduce supp/ even further, and this would contradict the choice of £. Therefore, 
e^T. □ 

Recall that given two subsets A and S of a metric space, the distance dist(j4, B) is defined by 
(27) dist(A, B) = sup(dist(a, B); a € A). 
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This function is not symmetric in general, and the Hausdorff distance is defined as 

distH(A, B) = max (dist(A, B), dist(B, A)) . 

Let r^i 3) J7 be another domain so that dili is given by dist(a;,9f2) = e <C 1. Let Es, si < s < S2 be a 
continuous family of compact oriented surfaces in Qi. Examples include also surfaces that are not closed in 
n but can be extended as closed ones in the larger domain fli. By a continuous family, we mean a family so 
that the Hausdorff distance distH(Ss, ^sg) tends to 0, as s sq, Vsq. We assume that each divides fti 
in two (open) parts: one, in the direction of the normal giving the orientation, that we denote by E"^'; and 
another one that contains dfli, that we denote by E™'. 
Let r C dfl be a relatively open subset of dil. Set 

(28) g ■.= {{t,x)- xeT,0<t<T{x)}, 

where r is a fixed continuous function on T. This corresponds to measurements taken at each x G F for the 
time interval < t < t{x). One special case is t{x) = T, for some T > 0; then Q = [0, T] x F. 




"inductive" step of the proof. The r.h.s. curve is Eg; the doted one to the left is at distance 
e from Es. 



Theorem 2.3. Let Qi, Q, and E^ he as above. Let u solve 
regularity conditions and the ellipticity condition 

(a) suppf C E|,f, 

(b) Es n O is strictly convex for any s, 

(c) ^,n{dn\r) ^% for any s. 

Assume that 



11) and have zero Cauchy data on Q. Let the 
13) be satisfied on K := (UE^) n fi. Assume that 



for any a; G E^. H fi, there is y Cz T so that T{y) > dist(a;, y). 



(29) 
Then 
(30) 

Proof. In this proof, we regard -F as a function supported in Q. Fix s g [si, S2], and assume that 



suppFnEs=0 Vs. 



(31) 

Then for any x e E., 
Proposition 



2.1 



suppF C E's"t. 

n n, t{x) > dist(x,y) for some y £ T. By the unique continuation statement of 
u has zero Cauchy data on E^ n f2 for \t\ ^ 1, and by assumption, this is also true on F. 
Let be a boundary normal coordinate to Eg positive in E'/"'. Let ipit, x) be as in (22 1, and e > be as 
in the proof of Theorem |2.2| The function ip — exp(/i0), 3> 1, is guaranteed to be pseudoconvex only for 
X in an 0{e) neighborhood of Eg n fj, and for \t\ — 0{y/e), if e ^ 1, by (b). We apply Lemma [2T] to the set 
D — E'"' n O to conclude by (b) and (c), that F = in some neighborhood of E^ n Ct, see Figure [2] 



Let now sq be the supremum of all s G [si,S2] for which (31) holds. The latter set is non-empty, by (a). 
By the continuity of s i-7> E^, (31 ) holds for s ~ sq. Indeed, assuming that (31 ) does not hold for s = sq, we 
can find < e <C 1, so that (31) does not hold for sq — e < s < sq. That contradicts the choice of sq to be 
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the least upper bound. On the other hand, by what we proved above, Sq cannot be an upper bound, unless 
s = 82- This completes the proof. □ 



Remark 2.1. The meaning of condition (29) is to guarantee that any point x where we want to prove 
F(x) = is reachable from F (from some point y) at a time not exceeding t(i/). In other words, there is a 
"signal" (a unit speed curve) issued from x that will reach the observation part F of d^l at a time while we 
are still making measurements there. By finite speed of propagation, it is a necessary condition, as well. 



Remark 2.2. A sufficient but an easier to formulate condition to replace (29) is 
(32) g=[0,T]xF with T > maxdist(S5 n 17, F), 



see (27). An even simpler sufficient condition is 

(33) g^[0,T]xr with T > dist(17,F). 

2.1. Examples. 

Example 2. Let fl C R^ be a bounded domain, and let g he a metric on Cl. Assume that dft is convex in 
the metric, and that there is xq G dil, so that all geodesies issued from dil, pointing into f2, exit after 
some fixed time. That property does not depend on the way we extend g outside Cl. All simple (see [H] for 
a definition) {^l,g) have this property, and all non-trapping convex ones have it, too. We will show that in 
this case we can cover by a foliation of smooth surfaces (curves, actually) Ss that are a small perturbation 
of the geodesies through xq. 



Figure 3. Example [2] 

Extend g in a small neighborhood fli of ft, and let yo ^ he close to xq so that the geodesies through j/o 
have the same property but in fli. Choose global coordinates in the latter as normal coordinates centered 
at yo- In those coordinates, the geodesies through yo £^re the lines, i.e., they solve x" = 0. Let at each x, J 
be the rotation operator in the tangent space given by Jv — {—V2,vi), where we used the standard index 
raising/lowering convention, and vector Jv is identified with the covector on the r.h.s. Given < (5 ^ 1, 
define the curves Ss as solutions of a;" = —SJx'. The parameter s measures the angle of the initial direction 
at yo with a fixed direction. Then are strictly convex when viewed from the side determined by the 
normal Jx' . In Figure [3) this is the upper side. We can always extend the curves Eg to closed ones with the 
extension being outside fl. Assume now that u has zero Cauchy data. Then we can apply Theorem |2.3| to 
conclude that F = when T is appropriately chosen. A possible choice of T is the diameter of fj, given by 
max(dist(a;, y); x,y Q Cl). To optimize T, we can consider a similar family, with 5 negative. In Figurejs] they 
are shown as dashed curves. Then the "interior" and the "exterior" are reversed. The two families converge 
to the set of the geodesies through yo, as \d\ ^ 0. Then the value for T enough to apply the theorem can 
be obtained by finding a geodesic 70 through yo so that maximum of the distances from 70 H Si to the upper 
and the lower side of dQ is maximized; then T is that value. 

Let he an ellipse, and let g he Euclidean. If yo is one of the vertices on the major (minor) axis, then 70 
is the major (minor) axis, and it is enough to take T to be a half of the length of the major (minor) axis. 
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Example 3. Let il C be as above. Assume that there exists a closed non-self-intersecting geodesic 70 of 
the metric g. Assume that the region between dfl and 70 can be foliated by a continuous family of strictly 
convex curves Sg. Then supp_F is contained inside 70, if T = dist(7o, 957). Our analysis does not allow us 
extend the equality F — inside. 




Figure 4. Example [3] 



Example 4. Let ft C be diffeomorphic to a disk, and let F be a relatively open connected part of dfl. Fix 
a metric g in some neighborhood of Ct. We do not need to assume that the whole dfl is convex but we will 
assume that there is a continuous family of geodesies, with endpoints outside 51, covering the region between 
F and the geodesic connecting the endpoints of F, see also Figure [5] In Figure [5] this is the unshaded region. 
The latter assumption is fulfilled if for example one of those points has the property that all geodesies issued 
from it, and pointed inside fl, leave the unshaded region after some fixed time. In particular, {^^g) being 
non-trapping suffices. Then we can perturb those geodesies to curves that are convex (bending to the shaded 
region), as in Example [2] to show that when T is chosen in an appropriate way, supp-F must be in the shaded 
region. 




Figure 5. Example |4] 

An higher dimensional analog of this example would be il diffeomorphic to a ball with F C dil diffeomor- 
phic to a disk on dft. Then = in the region covered by families of convex surfaces. For example, let 
be the ellipsoid Y^{x'f/a'f = 1, and let g be Euchdean. Let T = dnn{x^ > C} with < C < ai. Then 
F = Omnn{x^ > C}, and it is enough to choose T — ai~ C . That T may or may not be sharp, depending 
on all aj and C. 

3. A NON-LINEAR PROBLEM OF RECOVERY OF A SPEED WITH ONE MEASUREMENT. APPLICATIONS TO 

THERMOACOUSTICS 

In section [2j we showed that one can uniquely recover / when t varies over a symmetric interval [— T, T], 
and u\t=o is known. No knowledge of Ut\t=o is required. Assume now that t varies over the interval [0,T], 
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and a, u admit even extensions for t e [~T, T] of regularity as in the preceding section. In particular, this 
means that Ut\t=o — 0. In other words, u solves 

in (0,r) X 

(34) { u|t=o = m^l, 




m 



compare with ([2|. Then one has obvious corollaries of the results of the previous section that we are not 
going to formulate. One can interpret those results as a recovery of a source, given a. If a = 1, then one can 
differentiate the equation above w.r.t. t and to reduce the problem to recovery of an initial condition for the 



homogeneous wave equation, that is the classical thermoacoustic problem of recovering / given A/, see (35) 
below, with a known speed. 

3.1. The thermoacoustic model. Let u solve the problem ([T]) where c — c{x) > is smooth and T > 
is fixed. Note that the wave equation is solved in the whole R" now. 

Assume that / is supported in Cl, where C R" is some smooth bounded domain. Assume also that 
c — \ outside i7. This is not an essential assumption; it is enough c to be known and fixed outside f2. The 
measurements are modeled by the operator 

(35) A/ := u|[oj]xan- 

The problem is to reconstruct the unknown c and /, if possible. We will study now the case when / is 
known, and we want to reconstruct c. 

3.2. Uniqueness results for the linear problem. Let (c, /), (c, /) be two pairs, and let u, u be the 
corresponding solutions of ([T]). Then 

[dj -c^I\){u-u) = {c^-c^)Au in(0,T)xR", 

(36) { {u-u)\t=o = 0, 

dt{u-u)\t=o = 0. 



Then 

(A-A)/ = (S-u) 

We have 

(37) (5A)/:-(A-A)/ = 



[o,T]xaa" 



w\ 



[0,T]x9n' 



where w solves ^ with F and a as in ([s]). Then suppF C fi, and given the regularity of c and c, we 
also have F = on d^. The measurement ([37| however determines the Dirichlet data on [0,T] x dVl only. 



The Neumann data can be recovered from that however, see also [TBI sec. 7], and Lemma 3.1 below. We 
emphasize again that w solves the wave equation for x in the whole R", and this is what allows us to recover 
the Neumann data. 

We assume below that w solves the more general problem 

( {df-Ag)w = a{t,x)F{x) in(0,r)xR", 

(38) <^ w\t^a = 0, 

where g is a smooth Riemannian metric that is Euclidean outside fl. In some of the results below, we do 
not assume that F and a are given by 

We introduce the energy space associated with the wave equation below, to be able to deal both with 
metrics and variables speeds. Write the metric g as g = c~^go, where c{x) > 0. In the thermoacoustic case, 
one usually takes go to be Euclidean. Let Ag^ be the Laplace-Beltrami operator as above but related to g^. 
Modulo lower order terms, c^Ag^ and Ag coincide. Write P in the form 

(39) P = ~ A, A = c^Ag,, + lower order terms. 
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compare with Q. Notice first that c^^gf, is formally self-adjoint w.r.t. the measure c ^dVol. Given a 
domain U, and a function u{t^x), define the energy 

Ju 

where |Vj:u|o is the norm in the metric go, and d VoIq is the volume measure w.r.t. go as well. In particular, 
we define the space Ho(U) to be the completion of C^{U) under the Dirichlet norm 



(40) 



I/IIL 



|V^u|gdVolo(x) 



It is easy to see that Ho{U) C H^{U), if U is bounded with smooth boundary, therefore, Ho{U) is topolog- 
ically equivalent to Hq{U). li U = R", this is true for n > 3 only. By the finite speed of propagation, the 
solution with compactly supported Cauchy data always stays in even when n = 2. The energy norm for 
the Cauchy data [/i, /2], that we denote by || • is then defined by 



\[fj2Wn= I (|V,/i|2 + c-2|/2|2)dVolo(x). 
Ju 



This defines the energy space 



n{U)^HD{U)®L^{U). 
Here and below, £^(C/) = i^(C/; c^^dVolo). Note also that 

(41) II/IIL 

The wave equation then can be written down as the system 



(42) 



Ut = Au, A = 



/ 
A 



where u = [u, ut\ belongs to the energy space %. The operator A then extends naturally to a skew-self adjoint 
operator on T-L. We denote by U(i) the group exp(tA). In this paper, we will deal with either U = R" or 
U = 0,. In the latter case, the definition of H£i{U) refiects Dirichlet boundary conditions. 

We will define next the outgoing DN map. Given g e C^((0, oo) x dVl), let w solve the exterior mixed 
problem with c = 1: 



(43) 



Then we set 



(5? - A)t; 
'v\[o,T]^dn 
v\t=o 
dtv\t=o 



in (0, T) X R" \ O, 

9, 
0, 
0. 



Ng 



dw 
dv 



[0,T]x9n 



By 0, for g € i7(\,^([0,r] x dil), we have [w,wt] € C{[0,T); n{n)y, therefore. 



N : iJ(io)([0,T] X dn) C([0,T] x H^dfl)) 



is continuous, where the subscript (0) indicates the closed subspace of functions vanishing at t = 0. Note 
that the results in ^ require the domain to be bounded but by finite domain of dependence we can remove 
that restriction in our case. We also refer to [S] Proposition 2] for a sharp domain of dependence result for 
exterior problems. 

When F and a are given by (|3|, the next lemma follows directly from its version [161 sec. 7] for A/ by 
subtracting A/ and A/. 
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Lemma 3.1. Let w solve ^ with F supported in Vl, and let t f->- a(t, ■)F{-) G L'^{fl) be continuous. Assume 
that c ~ \ outside 51. Then for any T > 0, w\[Q x]xdn determines uniquely the normal derivative of w on 
[0, T] X dn as follows: 



(44) 



dw 
dv 



[o,T]xao 



T]xanj 



Proof. Let v be the solution of (431 with g — w|[o_T]xan- The latter is in i/^-'p^([0, T] x d^). Let w be the 
solution of ([2]). Then v — w solves the unit speed wave equation in [0, T] x R" \ f2 with zero Dirichlet data 
and zero initial data. Therefore, v — w va [0, T] x R" \ Vt. □ 



With this in mind, we have the following version of Theorem |2.2| in this context. Notice that in the two 
theorems below, we do not assume F and a to be given by ([s]). 

that dfl is strictly convex, w solves and the function a satisfies the elliptic 



Theorem 3.1. Assu 

condition [13] in the closure of the set (J^d) below. Let a and w admit even extensions satisfying (10). Let 
T > be such that a;" := dist(a;, dfl) is a well defined smooth function in 51 with non-zero differential; ana 
{x" = s} , < s < T are strictly convex surfaces. Then if w — on [0,T] x dQ, we also have 



F{x) ==0 for X en, dist(a:, dn) < T. 



(45) 

Proof. We first apply Lemma [3.1| to conclude that the normal derivative of w vanishes on [0,T] x 951, as 
well. The assumptions of the theorem imply that a(t, x) can be extended in an even way satisfying the 
assumptions of Theorem |2.2[ Then w admits an even extension, as a solution of the wave equation, and thus 
we apply Theorem 2.2 □ 



Since recovery of in (j2| from w|[o,t]x(50 is a linear problem, we also get uniqueness for that problem in 
the set ( 45 ) . We also get unique recovery of the speed c in the region in ( 45 ) . Those are in fact partial cases 



of Theorem 13.21 and Theorem 13.31 below. 

The local data result, in the spirit of Theorem 2.3 is not so straightforward because the recovery of the 
Neumann data is not so direct. 

Define the "cone" 



(46) 



{{t,x) e R+ X 951; t + disto(a;,i) < i}, 



where for a, b in R" \ 51, disto(a, h) is the infimum of the lengths of all smooth curves lying in R" \ 51 that 
connect a and b. 



Theorem 3.2. Let 51, 51i, and be as in Theorem 2.3 Let w solve fel) and let a, w have even extensions 



that satisfy (10) and {13) with K being the closure of the set in (47) below. Assume that w = on Q . Then 
(47) F = in {xenn (U,E,); 3y e F, Cdist(a:,y)..y C G}. 

Proof. Outside 51, w solves the wave equation Pw = with zero Cauchy data for t = and zero Dirichlet 



data on Q, see (28). This does not allow us immediately to conclude that the Neumann data vanishes there. 



too. On the other hand, by the finite domain of dependence result in [3J, dw/dv = near [t, y) £ G, if the 
"cone" Ct.y is contained in G- By Theorem 2.3 this implies w{x) = in a neighborhood of all x in the union 
of all so that dist{x,y) < t. This completes the proof. □ 



Remark 3.1. In case of observations on the whole boundary, i.e., when T = 951, with G = [0,r] x 951, (47) 
implies F = in the set dist(a;, 951) < T. In particular. 



(48) 



T > dist(51,951) 



is sufficient to conclude F = 0, see also (33). This in agreement with the results of the previous section since 



one can recover the Neumann data easily by Lemma 3.1 
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Example 5. Let g be Euclidean and let fl be strictly convex. Let F = {x" > a} H dCl with some fixed a, and 
let ila = {x" > a} n ri. Then ^la satisfies the foliation condition. For any x £ Hq, let y be the point on F 
with the same x' coordinates, where x' = (x^, . . . , x"~^). Then |a; — y| < disto(y, dil \ F) because even when 
X £ {x" = a} (then |x — yj = y" — a is maximized), the Euclidean distance from y to {x" — a} minimizes 
the distance from y to that plane with the constraint that we take it outside fl. Then the "cone" C\j:-y\^y is 



included in g, if we choose g = [0,T]xT with T = dist(F, dD, \ F), see (|27|). With that choice of g, under 
the assumptions of the theorem, we get = in O^. In other words, the "cone condition" in (47 1 is satisfied 
and therefore it is not restrictive in this case. 

By a perturbation argument, if g is close enough to the Euclidean metric, then F would vanish in a bit 
smaller set than fio- 

3.3. Uniqueness for the non-linear problem. We now go back to the problem of determining the sound 
speed c in ([I]) from A/ with / fixed and known. Clearly, some conditions on / are needed since when f — 0, 
for example, we get no information about c. 

Based on Theorem |3.1| or Theorem |3.2[ we can easily formulate versions for the non-linear problem. We 
will formulate a consequence of the latter theorem under conditions that guarantee that we can recover c in 
the whole fi. 



Theorem 3.3. Let c and c be two smooth positive speeds equal to 1 outside il. Let be as in Theorem 
and satisfy (a) and (b) there with F := (? — and with the strictly convexity assumption in (b) fulfilled 
w.r.t. the speed c. Let 



(49) A/ = A/ on [0,T] x dn, with T satisfying \3^j . 
Assume that for some compact K d fl, 

(50) supp(c -c)(lK, A/ ^ on K. 

Then c — c in US^. 

// in particular USs is dense in fJ, and T > dist(r2, d^), then c = c. 

Examples of the latter statement include geodesic balls with a fixed center under the assumption that 
they are all strictly convex. Then c — c everywhere except in the center, and by continuity, this is true in 
the center as well. 

Proof. Set w = u — u, where u and u solve ([T]) with the speeds c and c, respectively; and the same /. Then 
solves ([2]) with a and F as in ([3|. The condition (13) is then equivalent to Au|t=o 7^ 0. By ([!]), this is 



w 



equivalent to A/ 7^ 0. Then an application of Theorem 3.2 completes the proof. □ 



Remark 3.2. The condition A/ 7^ may look mysterious at first glance. The stability analysis below shows 
that it is needed for the linearization to be Fredholm. A simplified look at this condition is the following. 
Let us remove the need for c to be 1 outside Q. Then any harmonic function / is also a time independent 
solution M = / of the wave equation {df — c^(x)A)u — 0, regardless of c. Then A/ carries no information 
about c at all. If / is harmonic only on some open set U, then u = f (regardless of the speed) in the light 
cone with base U, and then a ~ there, see ([3|. Then the kernel of the linearized map would be C°° for 



y £U, as it follows from (53 1. That implies high instability for the linearization, at least. 



3.4. Stability. As a general principle, we have stability in Sobolev spaces if we can detect all singularities at 



g where wc make measurements, see, e.g., |13| . Under the assumptions of Theorem 3.3 that would require 
the following 

(51) V(x,^) G SK, the geodesic through (x,0 hits dVl for some t with \t\ < T. 

Here we used the fact that the problem extends in an even way w.r.t. t. We also identify vectors and covectors 
by the metric g. 



Notice that condition (51 1 is stronger than the uniqueness condition (48). The latter requires that from 



any x there is a signal (a unit speed curve) originating from x reaching 917 up to time T, i.e., dist(a:, dVL) < T. 



Condition (51 1 requires that from any x and any direction ^ the geodesic through it reaches 9f2 for time 
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\t\ < T. The same conditions appear in the analysis of the thermoacoustic problem of recovery of /, given c 
and A/, see [inilTI]. Here however, we assume the foliation condition as well. 

Let Rw be the trace of the first component it; of w := [w, Wt], defined on [0, T] x il, to [0, T] x dVl. 



Theorem 3.4. Let w solve ^ with w, a satisfying the regularity assumptions of Theorem 3.1 , and let F 
he supported in a compact K C ft, with F G L^{K), and let a(0, ■)\k ^ 0. Let Eg be as in Theorem 2.3 , and 
assume that UEs is dense in £7 (the foliation condition) . Let K and T satisfy \51^ (the stability condition). 
Then 

(52) WFWh'^i^K) <C\\wtt\\L^([o^T\xdn) 

with a constant C that remains uniform when the coefficient a stays in a fixed bounded set in C^([0, T]; C(f2)). 
Proof. We use the notation a{t) = a{t, •) below. Differentiate 

(53) w= / U(s) [0,a(t- s)i^] ds 
to get 

dtRw ^ RV{t)[0,a{0)F] + R [ V{s)[0,a' {t - s)F]ds. 

Jo 

Differentiate again and use the identity a'(0) = and the definition of A to get 

(54) dfRw = Aa(0)F + R [ U(s)[0, a"{t - s)F] ds = Aa(0)F + R [ U(t - s)[0, a"{s)F] ds. 

Jo Jo 

Let X e C°°([0,T]) be such that x = near T, and x = 1 on [0,To], where Tq < ^ is such that Q still 
holds with T replaced by Tq. Let B be the back-projection operator defined as follows. Let v solve 

{df-c^A)v = in(0,T)xl7, 
v\t=T = dtv\t=T = 0, 
v\[o,T]xdn = h. 

Then Bh := w|f=o- By [15j Theorem 3], Bx^ is a classical ^PDO of order with principal symbol 

where 7^^^ is the unit speed geodesic issued from (x,^), and ±t± > are the times need for it to hit dfl. 



We recall that we identify vectors and covectors by the metric g. Condition (51) guarantees that the symbol 
above is elliptic. Let Q be a parametric for it, i.e., QBx^ = Id + in fl, with K smoothing. For the 
purpose of this proof, we only need Kg g Vl/"^, that can be achieved with finite smoothness requirements on 
c. 



Apply QBx to ( 54 ) to get 



(55) QBxdtRw = {ld + Ko)a(0)F + QBx [ RV{t - s)[Q,a"{s)F]ds. 

Jo 

For any s G [Q,t], the function [0,a"{s)F] belongs to the energy space H and is supported in il. Then 
RU{t — s) (a more accurate notation would be i?U(- — s)) maps that function to a function that belongs to 
H^QjCR X 90), where the subscript (0) indicates a support disconnected from t ^ 0. This is explained in [15] 
in the context of thermoacoustic tomography, and the reason is that i?U(- — s) is an FIO of order with a 
canonical relation of graph type, see also |17) . The dependence on s is continuous, therefore, the integral in 
(54 1 belongs to that space, as well. By [3, i? : H^^s^CR. x 90) — Hjj{n) and is continuous. Therefore, the 
integral term in (54) is a compact operator of F in L^(0), mapping L^(0) into Hu{n,). 



Therefore, since a(0)|_R- 7^ 0, 
(56) \\F\\mK) < C\\d^Rw\\L-ao,T]xdn) + C\\K2Fh^^n), 

with K2 : L'^{K) — >• L^(0) compact. We used here the fact that B : L^Q,„p(R x 90) L'^{Vt) that also 
follows from ; or from the property of B to be an FIO of order with a canonical relation of graph type 
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[T5] . By [201 Proposition 5.3.1], estimate (56) implies a similar one, with a different C, and the last term 
missing. 

The statement about the uniformity of C does not follow directly from the last argument above because 
there is no control over C. Instead, we will perturb estimate (52 1. Notice first that by (54 1, the map 
C^([0,T]; C(il)) 9 a I— )■ Wtt\[o,T]y.dVi G L"^ is continuous. Then if a and a are two coefhcients and w, w are 
the corresponding solutions, we have 

\\F\\l-^(K) <C\\'Wtt\\L^([o,T]y.dn) 

<C'l|wtt||L2([0,T]xan) + C\\wtt - Wtt\\L^([O^T]y.dn) 

<C||wtt||L2([o,T]xao) + C5\\F\\l2(x), 

where 5^1 when a is close enough to a in C^([0, T]; C(f2)). We can therefore absorb the 5 term with the 
l.h.s. □ 

We are ready now to formulate the stability result for the thermoacoustic problem. 



Theorem 3.5. Let K C ft be a compact set. Let c, c € , f € H^^^ , k > n/2, be as in Theorem 3.3 and 
let the assumptions of that theorem be satisfied, except for (^^. Then 

(57) P-c||l2(k) <C||a2<5A/|U.([o,T]xaa), 

with C ~ C(Ci) uniform if c satisfies \\c\\(jk < Ci, 1/Ci < c, k > n/2. 

Proof. Apply Theorem 3.4 with a and F as in ([3|. We only need to prove the statement about the uniformity 



3.4 



of C. That requires to estimate ||u||c2([o,T]xOn) in terms of c, see ([3j) and Theorem 

It is straightforward to see that if c S C'=-\ then A'^f e H provided that f G H'^^x H''. Then A'=U(<)f 
is locally in H*^^^ x ij'^, therefore the first component of U(i)f is C''=+i-"/2 provided that fc + 1 — n/2 > 0. 
We have fc + 1 - n/2 > 2 when fc - 1 > n/2. Therefore, u £ C([0,T]; C^) when fc > n/2. The analysis of 
the rest of the second derivatives of u is similar. □ 
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